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Abstract 

It has recently been argued that the rate per unit volume of baryon number 
violation (topological transitions) in the hot, symmetric phase of electroweak 
theory is of the form rja^T^ in the weak-coupling limit, where t/ is a non- 
perturbative numerical coefficient. Over the past several years, there have 
been attempts to extract the rate of baryon number violation from real-time 
simulations of classical thermal field theory on a spatial lattice. Unfortu¬ 
nately, the coefficient rj will not be the same for classical lattice theories and 
the real quantum theory. However, by analyzing the appropriate effective 
theory on the lattice using the method of hard thermal loops, I show that the 
only obstruction to precisely relating the rates in the real and lattice theories 
is the fact that the long-distance physics on the lattice is not rotationally 
invariant. (This is unlike Euclidean-time measurements, where rotational in¬ 
variance is always recovered in the continuum limit.) I then propose how this 
violation of rotational invariance can be eliminated—and the real B viola¬ 
tion rate measured—by choosing an appropriate lattice Hamiltonian. I also 
propose a rough measure of the systematic error to be expected from using 
simpler, unimproved Hamiltonians. As a byproduct of my investigation, the 
plasma frequency and Debye mass are computed for classical thermal field 
theory on the lattice. 
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I. INTRODUCTION 


The violation of baryon number (B) in the hot, symmetric phaseQof electroweak theory 
plays a crucial role in scenarios for electroweak baryogenesis. The rate of B violation is tied, 
through the electroweak anomaly, to the the rate of topological transitions of the electroweak 
gauge fields. For the symmetric phase, this rate is not calculable by any perturbative 
method. In this paper, I address whether the topological transition rate can, in principle, 
be extracted from lattice simulations. The discussion will reveal some surprising features of 
real-time thermal field theory on the lattice. 

It has long been appreciated that, at finite temperature, topological transitions in real 
time are not directly related to topological transitions in Euclidean time P]. As a result, 
there is no apparent way to measure the real-time thermal transition rate in a standard. 
Euclidean-time, lattice simulation of quantum field theory. Several years ago, Ambjprn 
and Krasnitz cleverly implemented the observation [Q that a full simulation of quantum 
field theory is not actually required. Topological transitions occur in the symmetric phase 
through large configurations which are essentially classical. Indeed, all long-distance bosonic 
physics in a hot, ultrarelativistic plasma is effectively classical because of Bose enhancement. 
The number of quanta per mode in low-energy modes is given by the Bose distribution 

and, by the correspondence principle, this is the classical regime. Unlike quantum field 
theory, real-time simulations of classical field theory are tractable: just evolve the classical 
equations of motion. 

Classical thermal field theory has an ultraviolet catastrophe, historically famous in the 
context of black-body radiation. Because every mode has energy by the classical equipar- 


^ Here and throughout, I use the term “symmetric phase” loosely since, depending on the details of the 
Higgs sector, there many not be any sharp transition between the symmetric and “symmetry-broken” phases 
of the theory [^]J^. A sharp transition is in fact required for electroweak baryogenesis. 
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tition theorem, and because there are an inhnite number of modes per unit volume in con¬ 
tinuum held theory, the energy density is inhnite. In quantum held theory, in contrast, the 
ultraviolet contribution is cut oh at momenta of order T. Ambjprn and Krasnitz reasoned 
that the details of short-distance physics should not ahect the long-distance physics of topo¬ 
logical transitions. In their simulations, they put their classical system on a spatial lattice 
and then evolved it in continuous real time. The ultraviolet catastrophe was cut oh by the 
lattice spacing, which they progressively made smaller and smaller. 

Son, Yahe, and I have recently pointed out that short-distance ehects do not decouple 
as cleanly as Ambjprn and Krasnitz hoped. In particular, the short-distance modes cause 
damping of the long-distance dynamics, and this damping ahects the transition rate. We 
showed that damping reduces the transition rate by a factor of 0{a), where a is the weak 
hne structure constant. Because the damping is caused by short-distance physics, it is not 
universal, and so a comparison of theories with diherent short-distance physics becomes non¬ 
trivial. In the classical lattice theory used by Ambjprn and Krasnitz, for example, damping 
reduces the rate instead by a factor of 0{aaT) where a is the lattice spacing. As I shall 
discuss, there is an even more serious problem: because damping is dominated by short- 
distance physics, it knows about the anisotropies of the lattice. As we shall see explicitly, 
and as was hrst noted by Bodeker et al. 0, the effective long-distance physics of the classical 
lattice theory is not even rotationally invariant. This is in striking contrast to the familiar 
situation of Euclidean-time simulations, where rotational invariance is always recovered in 
the continuum limit. 

I shall assume that our analysis of damping in ref. 0 is correct and will eventually be 
borne out by numerical simulations at sufficiently weak coupling. In this paper, I focus on 
the natural follow-up question: Given that the real quantum theory and the classical lattice 
theory have different long-distance physics, is there any way to measure the real topological 
transition rate? The answer is yes—in principle—but an exact calculation requires a careful 
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choice of lattice action and the numerical extraction of a somewhat awkward limit .[| But 
I shall also propose a rough numerical measure of the suitability of generic lattice actions 
and argue that even results from simple, canonical actions should be in the right ballpark if 
properly interpreted. 

I shall not address at all the complicated problem of how one measures topological tran¬ 
sitions on the lattice in the first place, which has a long, difficult history for Euclidean-time 
quantum simulations^ and a shorter but still confusing one for real-time classical simula¬ 
tions |Hi|6|JTl|,[l^ ■ I simply focus on dynamics and assume the measurement problem will 
eventually be solved, using cooling or some other technique. 

Section 2 of this paper outlines my method for taking a measurement of the topological 
transition rate in a theory with one ultraviolet cut-off {e.g. a classical spatial lattice theory) 
and using it to predict the rate in the same theory with a different ultraviolet cut-off {e.g. 
the real, continuum, quantum field theory). The procedure will require, however, that the 
long-distance dynamics be rotationally invariant in both cases, and this poses a problem for 
the lattice that will eventually be dealt with. Section 2 rests on a very rough and schematic 
discussion of the effective long-distance dynamics, and I return in section 3 to do a better 
job of reviewing the correct long-distance theory. I review the derivation of “hard thermal 
loop” effective theory but add a small twist. The usual discussions in the literature are based 
on the assumption that hard {i.e. high-momentum) excitations in the plasma move at the 
speed of light. This is not true for lattice theories, and I show how the usual results easily 
generalize. In section 4, as a warm-up example of a rotational-invariant classical theory, I 
consider continuum classical thermal field theory with the ultraviolet regulated by higher- 
derivative interactions. In section 5, I turn to the canonical definition of the classical theory 
on a simple cubic spatial lattice. By explicit computation, the damping at long distances 

^ I shall briefly comment on the very different proposal of Hu and Muller in my conclusion. 

^ For a brief review, see sec. III.C.l of ref. |^. See also table I of ref. 
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is shown to be anisotropic and to have an interesting strncture of cusp and logarithmic 
singularities in its angular distribution. I discuss the origin of these singularities. Then I 
argue that measurements of the topological transition rate in simple lattice theories can still 
(if properly interpreted) be used to estimate the real transition rate, and I propose a rough 
measure of the systematic error arising from the anisotropy of the lattice. In section 6, I 
propose that the real transition rate could in principle be measured arbitrarily well from 
lattice simulations by implementing a lattice version of the higher-derivative continuum 
theory discussed in section 4. In section 7, I briefly discuss some possibilities of alternative 
lattice theories that may be more rotational-invariant than the canonical one, yet easier to 
implement than my proposal of section 6. Section 8 explores two quantities important to 
thermal held theory but slightly tangential to the main thrust of this paper: the plasma 
frequency and the Debye mass. I compute both for classical, thermal gauge theory on the 
lattice. The results provide one test of whether lattice simulations of a given size are indeed 
in the small coupling limit, and I comment on the application of this test to the simulations 
of ref. 1^. Section 9 offers my conclusions and a summary of what remains to be done. 

II. THE BASICS 

The origin of the basic scales associated with topological transitions is reviewed in the 
introduction of ref. [§], and I shall simply quote here the standard result: topological transi¬ 
tions proceed through non-perturbative gauge conhgurations of spatial size R l/g'^T. In 
ref. 0, we then showed that the time scale associated with such transitions is t 1/g^T 
and that the transition occurs through slowly varying magnetic conhgurations.^ (More tech¬ 
nically, the relevant, low-frequency gauge helds are transversally rather than longitudinally 
polarized.) My current endeavor can be summarized as follows. Consider the effective long¬ 
distance and long-time theory corresponding to frequencies and spatial momenta of order 


See also ref. [ p^ . 
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(a;, fc) ~ {g‘^T,g^T). What is the correspondence of those theories between (a) the real 
quantum theory, and (b) a cut-off classical theory? 

The important characteristics of the long-distance effective theory for topological transi¬ 
tions, as discussed in refs. [^|^, can be roughly introduced as follows: The behavior of the 
long-distance modes is analogous to a pendulum moving in hot molasses, where the molasses 
represents the short-distance modes, and non-perturbative physics such as the topological 
transition rate corresponds to large-amplitude fluctuations of the pendulum. A slightly 
non-linear pendulum is described by an equation of the form 

{d^ + k'^)A -|- (higher-order in A) = 0 . ( 2 . 1 ) 

A pendulum in (idealized) hot molasses is described instead by a Langevin equation of the 
form 


{dl -I- '^dt + k‘^)A + (higher-order in A) = ^{t). ( 2 . 2 ) 

There is a damping term with coefficient 7 due to the viscous interaction with the molasses. 
There is also a random force ^(t) due to the buffeting of the pendulum by random thermal 
fluctuations of the molasses. The two are related by the fluctuation-dissipation theorem. 
Ignoring the non-linearity of the pendulum, the random force in (p.2|) must have a white- 
noise spectrum with normalization 2'yT : 

= (2.3) 


The frequency response of such a pendulum is 


A{uj) = 


with a power spectrum 


{Aiu'r A{u;)) = 


uj'^ — i'^uj A k? ’ 
27 T 5{uj' — u) 


(2.4) 


(2.5) 


— a;2 — i'-yuj + 

which includes small-amplitude, high-frequency fluctuations and larger-amplitude, low- 
frequency ones. If the pendulum is strongly damped, 7 S> fc, then the characteristic inverse 
time for the largest amplitude fluctuations is 
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ijj ~ fc ^/7 <C k . 


( 2 . 6 ) 


In this case, the o;^ term is ignorable in (p.4|) , and the large-amplitude fluctuations of the 
system can be described by replacing (| 2 . 2 |) with the simpler equation 


{'^dt + k'^)A + (higher-order in A) = ^{t). 


(2.7) 


In the case of gauge theory, the damping coefficient^ 7 is determined from the imaginary 
part of the (retarded) thermal self-energy n(a;, k). The effective equation of the long-distance 
modes is, in frequency and momentum space. 


(—o;^ - 1 - n(a;, k) -|- k‘^)A -|- (higher-order in A) = ^{uj, k), (2.8) 


where k = |k|. The dominant contribution to 11 for a;,/c <C T is well-known |T^ and is 
generated by interactions with short-distance (momenta ~ T) degrees of freedom. In the 
frequency region of interest (a; -C /c -C T), and for transverse (he. magnetic) fluctuations, 
it is given in weak coupling by 


n(a;, k) —iu 


7rmj 


(2,9) 


where is the (leading order) Debye mass. For the case of pure gauge theory. 


ml = \CAg^T^ 


( 2 . 10 ) 


where Ca is the adjoint Casimir and Ca = N for SU(iV). So, for each mode k, the damping 
coefficient 7 of (p^) is just 


7 = 


Timl 


( 2 . 11 ) 


^ The reader should not confuse the damping coefficient 7 with the plasmon damping rate: they refer 
to different kinematic regimes. 7 describes the damping of nearly-static magnetic fields in the plasma, for 
which uj <^k and Imll is 0{g‘^T'^uj/k). In this case, Imll is dominated by hard contributions. The plasmon 
damping rate, on the other hand, describes the damping of propagating plasma waves, for which u > k 
and Imll is 0{g^T^) |l^. In this case, the dominant contribution to Imll depends on soft as well as hard 
physics. 
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For the spatial momentum scale of topological fluctuations, k = 0{g‘^T), ( p. 6 |) confirms 
that uj ~ g^T k and shows that the effective theory can be described by a simplihed 


Langevin equation of the form (|2.7| ). 

So far, I have been very schematic. Most signihcantly, I have blithely ignored all the 
terms in the above equations labeled “higher order in A.” Since I am ultimately interested in 
non-perturbative phenomena, I need to do a better job. I will review the full, unadulterated 
equations of motion of the effective theory in the next section. However, the schematic 
equations presented so far are adequate to explain the basic idea of how to compare theories 
with different ultraviolet behavior. 

Ignoring the higher-order terms altogether for now, consider two theories with different 
ultraviolet cut-offs—say, the real quantum theory and some cut-off classical theory. The 
self-energy H will be different in the two theories, and so the damping constant 7 will be 
different. The long-distance effective equations of motion are then also superhcially different: 


ilidt + k‘^)A — ^i{t) 


( 2 . 12 a) 


vs. 


(72^4 + k‘^)A — ^2{t) 


( 2 . 12 b) 


But, as I’ve written them, these two equations are trivially related by a rescaling of time: 

(2.13) 

72 

[That this transformation maps —>■ ^2 can be verihed from (|2.3|) .] That means that the 

topological transition rates T of the two theories are related by a simple rescaling of time: 

ri = —Ts. (2.14) 

7i 

The prescription for measuring the real topological rate is then: (a) to measure it in some 
UV cut-off classical theory, (b) to perturbatively compute the self-energy and hence the 
values of 7 in both the real theory and the classical theory, and then (c) to convert the 
measured rate by (|2.14|) . 






The fly in the ointment is that generically, in lattice theories, 7 is not rotationally 
invariant but depends on the direction of k with respect to the axis of the lattice. I shall 
demonstrate this later by explicit calculation. This means that there is no single rescaling 
of time that can relate the equations (| 2 . 12 | ) for all the long-distance modes of the theories. 
So there is another step required in the prescription: (step 0) And a classical lattice theory 
where the ultraviolet cut-off is as rotationally invariant as possible. I should emphasize that, 
unlike the more familiar case of Euclidean-time simulations, this last step is a theoretical 
necessity and not merely a numerical convenience. In Euclidean simulations, improving the 
rotational invariance of the action improves the rate of approach to the continuum limit as 
the lattice spacing is decreased, but (in principle) any lattice action will do as long as the 
lattice spacing is small enough. Here, that is not the case. 


The chink in the armor of rotation invariance can already be seen in the result (p.9|) for 
the self-energy in the real, rotational-invariant theory. In Euclidean time, consider effective 
interactions of the long-distance degrees of freedom that are induced by the short-distance 
degrees of freedom. Such interactions can be Taylor expanded in the small, long-distance 
momenta kx, ky, and kz- The cubic symmetry of the lattice is enough to guarantee that 
any interaction up to two derivatives k is in fact rotationally invariant as well. Interactions 
involving four derivatives need not be rotationally invariant, but such interactions are irrele¬ 
vant and decouple from the long-distance physics. The crucial assumption of this argument 
is the analyticity of the interactions at k = 0. It has long been known that real-time thermal 
interactions generated by short-distance physics fail this analyticity criterion, as can be seen 
explicitly in (|2.9| ). The standard mechanism by which rotational invariance is recovered in 
Euclidean time is therefore no longer operative. 

Before proceeding to flesh out the proposals made in this section, I should fix a point 
of nomenclature. Often in this paper I shall refer to the “real” rate Treai of topological 
transitions in continuum, quantum, pure non-Abelian gauge theory at high temperature 
and weak coupling. If one’s interest is in actual electroweak theory in the same limit, with 
its accompanying Higgs and fermion fields, it is easy to convert if one ignores the small 
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effects of the Weinberg mixing angle. The Debye mass in electroweak theory is |15 


= (I + bi) 


( 2 . 16 ) 


for a single donblet Higgs and rif = 3 families. By the conversion procednre (|2.14|) ontlined 
above and by (|2.11|) , the electroweak rate Tgw is related to the pnre SU(2) gange theory rate 


Treal” by 


p _ Tj-eal 

V4 2 


Uf 


(2.16) 


Eq. (p.l6|) and the analysis of this paper applies whenever the temperatnre is snfficiently 
high that the infrared dynamics of the Higgs is irrelevant at lengths of 0{l/g‘^T). This is 
the case either (1) far above the electroweak phase transition or crossover, or (2) in the 
symmetric phase at the transition in cases where there is a Ist-order transition and the 
transition is not exceedingly weak. Other cases could be handled by including the infrared 
dynamics of the Higgs into the discussion and into simulations. 


III. HARD THERMAL LOOPS 

A. Review of Basic Derivation 

It is now well-established how to write a set of non-local equations describing the domi¬ 
nant interactions of long-distance modes in hot, real-time, non-Abelian gauge theory when 
the random force term ^ discussed earlier is ignored.^ These interactions go by the name of 
“hard thermal loops,” since they are generated by integrating out the hard, short-distance 
modes of the theory. Though there is a great deal of literature on this subject, it is almost 


® The reader may wonder why numerical simulations don’t simply simulate these equations directly. The 
problem is their non-locality. For a discussion of the nightmare of trying to find a consistent, local description 
of the effective theory itself (which requires making precise the separation between soft and hard degrees 
of freedom), see ref. 0. I shall be using the hard-thermal-loop equations only as a theoretical tool for 
comparing the (leading-order) long-distance dynamics of more fundamental, local theories. 
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universally geared to hard particles that travel at the speed of light. In a lattice theory, in 
contrast, the hard particles have a more complicated dispersion relationship than uj = k. It 
will therefore be useful to briefly review the derivation of the hard thermal loop equations so 
that I can motivate the (simple) modihcations required for applying them to lattice theories. 

Hard thermal loops were originally derived through diagrammatic analysis by Braaten 
and Pisarski [l^. For my purpose, the quickest way to get at them is by an alternate method 
formalized by Blaizot and lancu |17|,[I^, which I shall present heuristically. (See also ref. [|T^.) 
The method is a generalization of the description of QED plasmas via Vlasov equations. The 
Vlasov equations are an effective description where soft photons are represented by gauge 
helds H^(t, x) in the usual way, but all the hard, charged particles are represented by a 
locally dehned flow density n(t, p, x), which is the density of particles at time t and position 
X with (hard) momentum p. The first Vlasov eq. follows from applying Liouville’s Theorem 
to n and using the equations of motion:]] 


0 = ^ n{t, p, x) = dtn + x ■ + p ■ dpU = + e(E + v x B) ■ dpU , 


(3.1) 


where v is the 3-velocity and where = (l,v). (Note that does not transform as a 
Lorentz 4-vector.) In the real world of ultrarelativistic plasmas, hard excitations move at 
the speed of light and v is just p = p/|p|. Here enters my simple generalization of the 
standard discussion in the literature: for more general underlying held theories, v should 
be identihed with the group velocity of the hard excitations, since it is the group velocity 
which describes the physical rate of motion of charge. If the dispersion relationship for the 
hard modes \s E = hip, then the group velocity is 


V — Vp Up . 


(3.2) 


(For those who hnd this style of argument overly heuristic, an explicit check of hnal results 
will be made in Appendix ^ using a diagrammatic approach.) 


^ For a discussion of why hard collision terms (he. terms non-linear in n) can be ignored in the application 
to topological transitions, see ref. Q. 
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The second Vlasov equation, which is for the gauge field, is simply 


= f = Y.'=i / Vn,, 


(3.3) 


where I have added a flavor index / to make explicit that there may be different types of 
particles with different charges in the plasma. I have also introduced the short-hand notation 

f _ f d^p 


( 2 ir)‘ 


(3.4) 


Now linearize the Vlasov equations in the gauge helds and in deviations 6n of n from its 
equilibrium distribution uq. Uq depends only on the energy hip of excitations and not on x 
or t. Using 


dno dno 

O'p^O = d'nii = V , 


dQ 


dQ 


(3.5) 


the linearized Vlasov equations are 


V ■ dduf -I- 6/ E ■ 


driQ 

dVt 


= 0 , 




(3.6a) 

(3.6b) 


(3.7) 


/ "P 

These equations can be solved formally for dn, yielding 

e is an inhnitesimal that has been inserted to pick out the retarded solutions to the equations. 
For QED with a single massless fermion, 

^ dno 




Q plJ-F ^ _4g2 / _ 

^ Jp dn {vd + e) 




(3.8) 


where the factor of 4 counts the two polarizations each of the fermion and anti-fermion. 

It is easy to guess how the forgoing generalizes to non-Abelian gauge theories. From 
( p.6b| ) it is clear that the 6n of interest must carry an adjoint color index. (Rather than a 
simple color-diagonal number operator, 6n now represents a correlation between the colors 
that the gauge bosons convert between.) But then the convective derivative n • cl in (|3.6aj) 
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must be generalized to a covariant convective derivative v ■ where D is understood to act 
in the adjoint representation. Similarly, the derivative in ( p.6b|) should be covariant. The 
resulting equation for the soft modes is 




drin 




V V- 


■E\ 


(3.9) 


J-p dVt {v ■ D + e) 

where has been replaced by C'a5 '^ and the factor of 2 counts the two helicity states of 
the hard gauge bosons. Readers desiring a much more detailed justihcation should refer to 


ref. [18 


For our current application to topological transitions, where frequencies are small com¬ 
pared to spatial momenta and the fluctuations are magnetic, we are only interested in the 


small frequency limit of (|3.9| ). Working in Aq = 0 gauge, this limit gives 

^ driQ 


1*9 j 


V V 




(3.10) 


where B = D x A. As long as the energy f2p respects parity, the integral over p would 
vanish if not for the e factor. Using 

1 „ 1 


V D -h e 


= P 


V D 


7r(5(v ■ fD), 


(3.11) 


where P stands for principal part, the effective equation for the frequencies and momenta 
of interest may be rewritten as 


dtA^ + (D X B)* = 0, 


(3.12) 


where 


7*^(k) = -27rCA5'^ 


driQ 

dVt 


5(v • k). 


(3.13) 


This is indeed of the form of the “non-linear system in molasses” ( |2.2| ) discussed schemati¬ 
cally in section |^, except that it is missing the random force term. 

The underlying diagrammatic origin of the result (|3.13|) for the damping factor can be 
roughly understood by considering the process of £g. |I|, which shows a soft excitation being 
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absorbed by the small-angle scattering of a hard excitation. The factor g n* in (|3.13|) comes 
from the interaction vertex, and the second factor gv^ is because the amplitude of £g. 
should be squared to get a rate. The (5(v ■ k) is just the uj k p limit of the energy- 
conserving delta function 5(r2p+k — hip — a;) associated with this diagram. Finally, dn/dfl 
for a Bose or Fermi distribution is just —(3 times the product no(l ± uq) of the initial-state 
probability distribution and the hnal-state Bose enhancement or Fermi blocking factor. More 
details are given in Appendix 


p hard 



p+k hard 


FIG. 1. The absorption of a soft excitation (wavy line) by the small-angle scattering of a hard excitation 


(straight line). 


B. The random force term 


The hnal result (|3.12|) of the previous section is in fact sufficient for the comparison of 
different theories, but, for the sake of completeness, I will review the random force term as 
well. The Vlasov equations of the previous section describe the relaxation or propagation 
of deviations from equilibrium; they do not, as formulated above, describe the fluctuations 
that occur in equilibrium itself. However, fluctuations are related to dissipation through the 
fluctuation-dissipation theorem. Introduce a random force term by writing 


7 *^'(iD) dtA> + (D X B)* = C • 


(3.14) 


If we momentarily ignore the non-linearity of (|3.14|), which arises from the non-Abelian 


nature of the theory, then the characteristics of the random force can be taken from (|2.3|): 


f“(t,x) T5®(x - x') 5{t - t '), 


(3.15) 
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where a and b are adjoint color indices. Hnet and Son have shown how to modify this 
eqnation to include the non-linearity and restore gauge invariance. The result is to simply 
replace 5“^ above by a straight, adjoint-charge Wilson line connecting x to x': 


r(t, x) x')) = 2 X, x') 7«T i<=»(x - x') i(( - «'). 


(3.16) 


where 


U{t; X, x') = V exp — / dz ■ A“(t, z)T“ ) , 


(3.17) 


z runs on a straight path from x to x', V denotes path ordering, and the are the (real, 
anti-Hermitian) adjoint representation generators. The basic equations for the effective 
theory are then (|3.14|) and (|3.16|) . 


C. Application to rotational-invariant theories 


The damping coefficient 7 of (^.131) can be rewritten as 

7rM*^'(k) 


7*^(k) = 


4|k| ’ 


where 


M*^(k) = -SC^g^ 


duQ 

dVt 


v'^v^ 5(v • k) 


depends only on the direction k. Note that M is transverse: 


(3.18) 


(3.19) 


i:‘Af9(k) = 0. 


(3.20) 


In a theory that is fundamentally rotationally invariant, it must then have the form 

M^^(k) = (W - k^y)m^, 

where is a number independent of k, obtained by averaging trM over k:R 


(3.21) 


I am using a notation and normalization here that makes m the analog of the in (2.9). However, m 


is not actually the Debye mass in theories where |v| 7 ^ 1. See section VIII. 
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(3.22) 


77^2 ^ -2C'a5'^ 


driQ 

dVt 


By examining the effective equations (|3.14|) and (|3.16|) , it is then easy to check that the time 
rescaling suggested in section || will indeed relate any two rotational-invariant theories: 


mi 

ti = — 5 -12 , 

mi 


Bi = —2 ^2 • 


(3.23) 


If flp is monotonically increasing, can be rewritten in the simpler form 

2 2 C'a ^2 00 

m = - — dppno{ilp). 

71^ Jo 


(3.24) 


The result ( p.lO|) for the real quantum theory can now be reproduced by using the massless 
Bose distribution (e^^ — 1)“^ for no(f2p). 


IV. AN EXAMPLE: CONTINUUM THEORY WITH HIGHER DERIVATIVES 


Before I move on to lattice theories and the difficulties with rotational invariance, it is 
useful to first make the preceding discussion concrete with a simple, rotational-invariant 
example of a classical thermal field theory. Consider a classical continuum description of 
non-Abelian gauge theory where I regulate the ultraviolet by higher-derivative terms.^ The 
Hamiltonian density is 


H = tT 


+ B^ f 



5* 


(4.1) 


where / has an expansion 


/ {z) — 1 -|- CLiZ -|- Cl2z'^ 


(4.2) 


® The term “higher derivative theory” often refers to Lorentz-invariant theories, where the action must con¬ 
tain higher time derivatives as well as higher spatial ones. Such theories create annoying issues about initial 
conditions and negative-norm states. The theory here, however, does not involve higher time derivatives 
and is free from such complications. 
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and E = —dtA is (minus) the conjugate momentum to A. At momenta low compared to 
A, this is just normal gauge theory. However, provided f{z) grows sufficiently rapidly at 
large 2 ;, the higher-derivative interactions will cut off the ultraviolet catastrophe of classical 
thermal held theory at a momentum scale of order A. 

To compare this cut-off classical theory to the real quantum theory, I will compute the 


damping coefficient of (|3.24|) . For any classical theory, the equilibrium distribution uo(hl) 
is simply 




(4.3) 


One quick way to see this is to start with the Bose distribution of the quantum theory and, 
for the hrst and last time in this paper, put in the explicit factor of h: 


no(0) = 




T 


as 




(4.4) 


This just recapitulates the discussion (|1.1|) of the introduction. The remaining element we 
need is the frequency Op as a function of momentum p for perturbative excitations. For the 
theory (|0|), it is 




A2 


(4.5) 


The coefficient of (|3.24|) is then (for monotonically increasing /) 


= 


2CAg^TA ds 


TT^ 


fis^) 


(4.6) 


This converges provided f{z) grows faster than z for large z. In other words, a four-derivative 
interaction is not enough to cut off the ultraviolet, but six derivatives do the job. By (|3.23|) 
and ( p.lO|) , the topological transition rate Fh.d. in this classical, higher-derivative theory is 
related (in weak coupling) to the rate Fj-eai in the real (pure gauge) quantum theory by 


kreal — T 


6 A ds 


h.d. 






(4.7) 


It will be useful later on to have the specihc result for / of the form f{z) = 1 -|- 
corresponding to an isolated six-derivative cut-off. In this case. 
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(4.8) 



where r(^) is the usual Euler Gamma function. 


V. A SIMPLE CUBIC SPATIAL LATTICE 


A. Results for 7 *^ 


We now have almost everything we need to compute the damping coefficient 7 ®-^, or 
equivalently for lattice theories; we just need the dispersion relationship Gp. Let’s 
start by studying the simple, classical, Kogut-Susskind Hamiltonian used by Ambprn and 
Krasnitz for simulations on a simple cubic (SC), spatial lattice. The tree-level dispersion 
relationship for propagating excitations on the lattice is 



(5.1) 


where a is the lattice spacing. For thermal held theory on the lattice, the parameters g^, 
a, and T always appear in the combination g‘^aT for dimensionless quantities.|^ I shall 
henceforth work in lattice units by setting a = T = 1. The conversion back to physical 





(5.2) 


Plugging into (|3.13|) with the classical distribution (|4.3|) then yields[3 


This is because and a can be scaled out of the Hamiltonian H so that the distribution exp(—/3R) 
becomes exp{—H/g‘^aT). 

Lattice papers such as 0 typically refer to the quantity /3l = ‘2N/gf^^ for SU(A) rather than gf^^.. 


A general formula for the hard contribution to H^ in a simple cubic lattice theory was first written 


down for QED coupled to massless scalars in eq. (Bll) of ref. 0. My (^) is that equation’s w <C fc limit, 


once the overall normalization in ref. [Q is changed to be appropriate for pure, non-Abelian gauge theory. 
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(5.3a) 


7*^'(k) = 


7rM*^'(k) 


M*^(k) = SC Kg 


4|k| ’ 

2 /■+’^ d^p sin Pi sinpj 




I-. (27r)3 

where the integrals are over the Brillonin zone \pi\ < vr. 


5(k ■ sin p), 


(5.3b) 



(b) 


(c) 


FIG. 2. The trace 7 ” of the damping coefficient as a function of direction, normalized by its angular 
average, for (a) a rotationally invariant theory, (b) a Kogut-Susskind Hamiltonian on an simple cubic lattice, 
and (c) for a face-centered cubic lattice. The spikes in (b) and (c) represent logarithmic divergences. 


The delta fnnction can be nsed to eliminate one of the three integrations, and I have done 
the remaining two nnmerically. (The detailed expression is given in Appendix |^) Fig. |b 
shows the resnlting angnlar distribntion of the trace 7 “ normalized by its angnlar average 
( 7 **). Not only is the resnlt not rotationally invariant, bnt it has an interesting strnctnre of 
singnlarities as well. There are cusps and, in certain special directions, 7 *-^ is inhnite! As we 
shall see, these singularities are mild and integrable. Also, they are true singularities only 
in the weak-coupling limit: they are rounded off by higher-order corrections but become 
sharper and sharper as p — 0 . 


B. The origin of singularities 


For any hxed direction k of the soft momentum, 7 *-^ as given by (p.l3|) can be rewritten 
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7*-’ (k) = -27iCAg‘^ 


(5.4) 


r 

Is dVL 




V V' 


|Vp(v 


k)r 


where the integral is over the surface S' in p space on which v(p) is orthogonal to k. If there 
is a point on this surface where the vector Vp(v • k) vanishes, there will be a singularity 
in 7. If there is a entire line of points where the vector vanishes, the singularity will be 
more severe. This is analogous to Van Howe singularities in the density of states.f^ Naive 
counting of constraints suggests the existence of singularities is to be expected: 5 degrees of 
freedom in p and k versus 4 constraints in v ■ k = 0 and Vp(v ■ k) = 0. 

To summarize, singularities in 7 will be caused by momenta p such that there is a null 
vector k, orthogonal to v = Vflp, of the matrix Vv = VVflp. In the case of the dispersion 
relationship (|5.1|) , 




(5®-^ cos p® — 


(6.6) 


At p = { 71 / 2 , 71 / 2 ^ 7 ^/ 2 ) and its lattice reflections, this matrix is proportional to u®u-^ and so 
every direction of k orthogonal to v cx (1,1,1) is associated with a singularity of 7. This is 
the origin of the cusp singularities in £g. |p3, which are along the plane kx + ky + kz = ^ and 
its three lattice reflections. 

For p = (7r/2, 7t/2,Pz) and general fhe matrix (|5.5|) has (1,-1,0) as a null vector, which 
is orthogonal to v. Therefore, k in the (1,-1,0) direction or its lattice reflections is associated 
with an entire line of singularities in p space. This is the origin of the inhnities in these k 
directions in £g. ^3. (See section |VII B| for a more general discussion.) For these values of k, 
Vp(v ■ k) vanishes linearly as one approaches the p-space line of singularities. The surface 
integral in (|5.4|) is therefore only logarithmically singular. 

Physically, there are not, in fact, inhnite spikes in the angular distribution of Imn(a;, k) 
in the lattice theory. Recall that {oJ,k) ~ (p^jP^) for the long-distance physics of interest. 


The density of states is p{E) = JpS(E — lip) and has singularities (typically cusps) when v = VPp 
vanishes for some p. 
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The infinite size of the spikes occurs only in the formal limit a; -C fc -C 1 that I applied 
to study the theory in weak coupling. The taming of the singularity can be seen from the 
Vlasov equations by returning to the hard-thermal loop equation ( |3.9| ) before extracting the 
small 10 limit. The imaginary part of the self-energy, for spatial polarizations, is 


Imn*^ = 


driQ 

dVt 


v^v^5{v ■ K ), 


(5.6) 


where K is the four-vector (cn, k). The relevant difference from the earlier derivation of the 
damping constant of (|3.13|) is simply that 5(v ■ k) has been replaced by 5(v • k — a;). 
This difference is sufficient to cut off the logarithmic infinities in from order g^k~^ In oo 
to order g^k~^ ln(fc/a;) ~ g^k' -Mn(l/^2)_ 

It is beyond my present purpose to explore the cut off of the logarithmic singularities in 
detail, but I should mention that the previous paragraph is misleadingly incomplete. Based 
on the Vlasov equations. I’ve shown how the logarithms are cut off by effects sub-leading 
in the limit a; -C fc <C 1. At sub-leading order, however, the Vlasov equations themselves 
are inadequate. One can see this from a more fundamental, diagrammatic analysis of the 
self-energy, discussed in Appendix The 5{y ■ K) discussed above is nothing more than 
the oo,k ^ p limit of the energy-conserving delta function — Op — a;) associated 

with fig. D and discussed at the end of section [HI A| . For (a;, k) of order {g"^, g^) in coupling, 
however, the correct sub-leading version is 


^(i^p+k — cUp — cn) 5(v • k — a; -f I k ■ Vv ■ k). (5.7) 

The k ■ Vv ■ k term is the same order as the oo term. Physically, it represents dispersion 
of the wave packets of the hard excitations—an effect not built into the Vlasov equations’ 
implicit treatment of hard excitations as a collection of classical particles. 

In any case, the divergences in the small-coupling limit are only logarithmic and so 
are mild enough not to cause any real problem with the effective theory specified by (|5.3bD . 
Topological transitions occur through configurations with spatial size 1/g‘^T and not through 
pure, infinite-extent, plane waves with momenta k precisely in one of the divergent directions. 
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Since logarithmic singularities are integrable in k, they will give only a hnite damping effect 
to the evolution of any spatially localized configuration. 


C. Estimating the real transition rate 


If there were no infinite spikes, fig. would look like at least a rough approximation to 
the rotational-invariant sphere of fig. ^ja, and a procedure for estimating the real, quantum 
transition rate Freai suggests itself; replace 7 by its appropriate angular average 7 , and use 
this value to relate the lattice measurement of F to F^eai as in (|2.14|) and section |III Q . What 
would be the expected systematic error of this procedure? Suppose 7 max(^) and 7 min(^) 
were the maximum and the minimum of the angular distribution of 7 (k). (I shan’t worry 
about the indices i,j on 7 ®-^ for the moment.) Since the effect of damping decreases the 
transition rate, it seems reasonable that the measured F is bigger than it would have been 
if 7 (k) had been isotropically equal to 7 max ia all directions but smaller than if 7 (k) had 
been isotropically equal to 7 min in all directions. So the ratio of the extreme values 7 min 
and 7 max to 7 would give a conservative estimate of the relative systematic error due to the 
anisotropy of the damping. 

I shall discuss later whether one might find lattice theories for which 7 (k) indeed has no 
infinite spikes. For the moment, focus on the case at hand. Since the spikes are integrable, we 
can still replace 7 by its average to obtain an estimate of F^eai from the lattice measurement 
of F. 7 is the damping felt by isotropic gauge configurations. The problem that remains is 
to find a plausible estimate of the systematic error, so that there is a measure for deciding 
whether one lattice theory with spikes is better than another. One obvious candidate is the 


root-mean-square deviation given by cr^ = (7 — 7 )^, where the over line indicates angular 
averaging. 

I have suppressed spatial indices on 7 above and will now be more concrete. By com¬ 
parison with the rotational-invariant case, where 7 *-^(k) has the form 


7®^(k) = (W-Ffc^)7(fc) 


(6.8) 
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due to its transversality, define 


f^'(k) = {5^^ -ley)^{k), 
7(*:) ^ Hy(k)>t, 


(6.9) 

(5.10) 


for the lattice case, where (■ ■ ■)j^ denotes averaging over the direction of k. The variation is 


a 


= i{7«(k)y(k))fc-|7(l;)p. 


(5.11) 


Using the general result ( p.l3| ) for 7 *-^ in any theory, we have 

_ ttCaS'^ /■ dno I I 

^w=-—ypd!Ti''i 


(5.12) 


and 


a 

^2 


driQ driQ 4(v ■ v') 


/\2 


pjp' dkl dVt' 7 r|v X v'l 

' [ ^ IviV 
Jp dVt I 


- 1 


(5.13) 


For the Kogut-Susskind Hamiltonian on a simple cubic lattice, numerical integration 
gives (j /7 = 0.31 and 


^{k) = 0.2687- 


k 


(5.14) 


Using (|3.23|) and ( p.ll| ), one can then estimate the real transition rate (in the weak coupling 
limit) from a lattice measurement as follows. The B violation rate T is expected to scale as 
in weak coupling |P]. In lattice units, measure the proportionality constant 

Tlat 


Then the estimate is 


rji = lim . 

Siat^O 


12 X 0.2687 

r..ai ~- Vi9 • 


(5.15) 


TT 


(5.16) 


The fact that the rough indicator cr/y of the systematic error due to using the lattice is 
about 30% suggests that the estimate ( b.l 6 |) should at least be in the right ballpark of the 
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true answer. (Also, the minimum value of tr 7 (k) in fig. Ipj is just 30% below the average.) 
Unfortunately, rji has not yet been measured,|3 and even the weak-coupling scaling law 
r ~ a® has yet to be verihed on the lattice. 


VI. ACHIEVING THE ROTATIONAL-INVARIANT LIMIT 


To eliminate the fundamental systematic uncertainty caused by the anisotropy of the 
lattice, one needs to hnd an effectively rotational-invariant lattice theory. (I re-emphasize 
that the search for an alternative theory is a theoretical necessity: just taking the small 
lattice spacing limit is by itself inadequate.) As it turns out, an effectively rotational- 
invariant theory is, in principle, easy to achieve. The idea is simply to start with a rotational- 
invariant continuum theory with a higher derivative cut-off, such as discussed in section B 
and then put it on a lattice that is hne compared to the continuum cut-off scale A“^. A 
numerical extraction of T would then require a careful double limit; 


lim lim Tiat 

aA^O giat^O 


( 6 . 1 ) 


where the order of limits is crucial. 

To be a little more specihc about how the cut-off A might be achieved on a spatial lattice, 
consider the tree-level dispersion relationship for propagating gluons with momentum small 
compared to the lattice spacing. For a given lattice Hamiltonian (with cubic symmetry), it 
will have an expansion in momentum of the form 


sSJp = ()|pP+ c,|p|* + C 2 (pJ+pJ + p*) +0(p“) 


( 6 . 2 ) 


in lattice units, where the coefficient s is just a reminder that the units of time can be 
normalized arbitrarily. By including next-to-nearest-neighbor couplings, and perhaps next- 
to-next-to-nearest neighbor couplings, and so forth, the parameters of the interactions can be 


rji should not be confused with the coefficient k presented in ref. jj], which was an attempt to extract 
the coefficient of a presumed F ~ scaling law for weak coupling—a scaling law which ignores the effects 
of damping 
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tuned to make any finite set of the coefficients 6 , ci, C 2 , • • • of the expansion take on whatever 
values desired. In particular, there is some local Hamiltonian for which (1) 6 is a small 
number, which I’ll call e^, (2) ci is 0(1), (3) C 2 = 0, and (4) all other coefficients are < 0(1): 

sHp = e^lpp + CiIpI"^ + 0(p®). (6.3) 

Choosing units of time where s = e^, 

^p = |pr + ^|p|^ + o(^/) . (6.4) 

This describes a theory with a rotational-invariant four-derivative “cut off” at momentum 
A = ea~^, plus anisotropic terms that are suppressed by at that scale. 

Unfortunately, the four-derivative “cut-off” in (|6.4|) is inadequate. As discussed in sec¬ 
tion |I^, a six-derivative interaction is required to cut off the UV contribution to damping. 
So, even for small e, a generic theory with dispersion relationship (|6.3|) will still produce sig- 
nihcant anisotropy.[^ We can £x the cut-off by going to six derivatives, requiring improved 
Hamiltonians where (1) 6 is a small number, which I’ll now call e^, (2) the fourth-order 
coefficients ci and C 2 vanish, (3) the coefficient of |p|® is 0(1), (4) the sixth-order coefficients 
that break rotational invariance vanish, and (5) other coefficients are < 0(1): 

sHp = di|p|®0(p®), (6.5) 

Choosing units of time where s = e^, 

= |pr +^|p|® + 0 . (6.6) 

This describes the desired rotational-invariant six-derivative cut-off at momentum A = ea~^ 
plus anisotropies that are suppressed by at that scale. The relative size of the anisotropies 
in the damping coefficient 7 will therefore be O(e^) and disappear in the limit e —> 0. 


The relative contribution of p ^ 1 to the damping rate is actually 0(1/In e) and so, technically, is small 
in the limit that In e was very large. But this is certainly an impractical limit for any numerical simulation, 
keeping in mind that, as discussed below, the e ^ 0 and t/iat ^ 0 limits do not commute. 
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I will not attempt in the present work to explicitly construct the improved lattice Hamil¬ 
tonian which gives this behavior. A very nice discussion of improving lattice Hamilto¬ 
nians is given by Moore in ref. |^, where he constructs a classical Yang-Mills Hamiltonian 
with vanishing terms. Constructing a Hamiltonian that satishes the other constraints for 
is presumably straightforward in principle but tedious in practice. 

In everything above, I have referred to conditions on the tree-level dispersion relationship. 
The real dispersion relationship for hard excitations {p ~ A) will be modihed by perturbative 
corrections, whose strength is parametrized by giat/Aa = (^lat/e. This is the reason that 
taking the gi^t —>■ 0 limit before the Aa ^ 0 limit would be crucial to extracting results from 
numerical simulations. 

For an improved Hamiltonian giving the desired low-momentum dispersion relation ( |6.5|) , 
the conversion of lattice measurements to the real transition rate would be achieved by hrst 
measuring the coefficient 

£ hiat 


P 2 = lim lim 


10 




(6.7) 


and then taking 


p _ 

real 






( 6 , 8 ) 


2 

which follows from (14. 8 |) . gi^^ above continues to mean g‘^aT where g^ is the continuum 
coupling extracted from the Hamiltonian . 


VII. OTHER LATTICE THEORIES 

A. The FCC lattice 

The aA —> 0 limit discussed in the last section may be difficult to extract numerically. 
So it’s worthwhile to consider whether there are relatively simple lattice theories where the 
damping coefficient 7 , though anisotropic, is closer to rotational invariance than for the 
simple lattice Hamiltonian discussed in section ^ I will briefly discuss here one possibility 
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for improvement: replacing the simple cubic lattice by a face-centered cubic (FCC) lattice.|^ 
The FCC lattice can plausibly lead to more rotational-invariant results because sites have 
more nearest neighbors than for the simple cubic lattice (eight instead of six). 


The formulation of lattice QED on a spatial FCC lattice is discussed by Gosar in ref. [23 


The fundamental plaquettes of an FCC lattice are equilateral triangles, and Gosar studies 
a simple Kogut-Susskind Hamiltonian dehned on those plaquettes. All I need here is the 
dispersion relation, which will be the same for a non-Abelian gauge theory as for photons in 
QED and is given by Gosar. For a simple cubic lattice, there are three orientations of links 
(h, y, and z) and three branches of excitations: two, degenerate, propagating polarizations 
with energy (|5.1|) ; and one non-propagating polarization with Dp = 0. For an FGG lattice 
there are six orientations of links and so six branches of excitations [^: 


(a) al = o, 


(7.1) 

(b) Ql = 32, 


(7.2) 

c,d) nl = 32 

1 - ^1 - ep/16 

(7.3) 

(e,f) nl = 32 

l + C-'p/16 

(7.4) 


where Cp is the energy a massless scalar held with nearest-neighbor coupling would have on 
an FGG lattice. 


ep = 4 


„ Px Py Py Pz Pz Px 

6 — COS — COS — — COS — COS-COS — cos — 


(7.6) 


2 2 2 2 2 2 

I have used lattice units where the side of a unit cell is a = In the long wavelength 
limit, branches (c,d) correspond to the two polarizations of the continuum case D = |p|. 


In four-dimensional Euclidean gauge theories, people have tried different lattice types for the merely 
practical reason of trying to speed approach the continuum limit. Lattices used include the F 4 lattice HI, 
which is the four-dimensional analog of an FCC lattice, and the body-centered hypercubic (BCH) lattice 

H- 

Ref. instead uses units where the distance between nearest neighbors is 1. 
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The non-propagating branches (a) and (b) do not contribute to the damping coefficient 7 
because their group velocity is zero. 

The shape of the hrst Brillouin zone of an FCC lattice is slightly complicated (a “trun¬ 
cated octahedron”), but all that is needed here is that the cube \pi\ < 2ti in momentum 
space covers the unique momentum states exactly twice. So, when integrating over states, I 
take 


1 d?p 

2^i-2,r (27r)3 ■ 

The sum is a notational reminder to sum over the two dispersion relationships 



(7.6) 


f!p = 32 


l±\/l-ep/16 


(7.7) 


The corresponding group velocities are 


= ± 


■ Px Py . Pz 
sm — cos — S- cos — 

2 V 2 2 


(7.8) 


- ep/16 

and its permutations. 

Numerical evaluation of the general formula ( 3.13 ) for 7 (see Appendix P) gives the 
angular distribution of tr 7 shown in £g. |^c. The position of the logarithmic spikes is 


explained in the next section. Numerical integration of ( 5.12 ) and ( 5.13 ) gives ^{k) = 
0.501(1) /k and cr-y/'^ = 0.26(1). The 26% result for the angular deviation of damping 
on the FCC lattice is indeed better than the 31% result for the SC lattice, but not by a lot. 
The minimum of tr 7 (k) in hg. |c is 20 % below the average. 


B. Avoiding logarithmic spikes 


The simplest actions on both the SC and FCC lattices have logarithmic spikes in the an¬ 
gular distribution of 7 , shown in figs. |b and c. In the proposal of section ^ for reaching the 
rotational-invariant limit, it doesn’t matter whether the hardest modes p ~ 1 of the theory 
give rise to such spikes, because the effect of the hardest modes, and hence the strength of 
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the spikes, will vanish as e —0. For the purposes of surveying simpler lattice theories that 
are only approximately rotationally invariant, however, it’s interesting to consider whether 
it’s possible to avoid the spikes altogether. So I will take a moment to explain the generic 
property of the simple SC and FCC actions that causes the spikes and then discuss whether 
it can be avoided. 


Recall the conditions k ■ Vfl = 0 and k • V Vfl = 0 of section [VB| for generating some sort 
of singularity. Consider the energy flp of excitations for small but hxed pz in a lattice theory 
with cubic symmetry. Fig. |a shows a qualitative contour plot of energy flp vs. px and Py 
for some such lattice theory [e.g. the simple cubic theory of (|5.1| )]. There is a minimum at 
Px = Py = which there must be if we are plotting the branch of flp that approaches the 
continuum limit flp = |p| for small p. Suppose that there is also a maximum somewhere 
along the line Px = Py, as shown in £g. ^a. Then one can argue as follows that there must 
be a singularity associated with the direction k = (1, —1,0) and with p somewhere on the 
line Px = Py {pz still hxed). By rehection symmetry through the plane Px = Py, the group 
velocity Vfl on the line must be perpendicular to k, and the curvature VVfl on the line 
must have the block-diagonal form 


(7.9) 



(A 

B 

o\ 

VVfi = 

B 

A 

0 


lo 

0 

cj 


in the basis {pz, Px+Py, k) • Near the minimum in the hxed pz plane, C must be negative. Near 
the maximum that by assumption lies on the line px = Py, it must be positive. Therefore, 
there is a point inbetween where C = 0 and so k ■ VVff = 0. This point then satishes both 
conditions for generating a singularity associated with the direction k = (1, —1, 0). 


29 






(a) 


(b) 


FIG. 3. A qualitative contour plot of flp vs. px and py for small, fixed pz in two different lattice theories. 
Black corresponds to low values of flp and white to high ones. 


But now note that this argument did not depend on the value of As long as the 
assumption that both a minimum and a maximum in {px,Py) simultaneously lie on the line 
Px = Py holds for some range of Pz, then there will be an entire curve of points in p space 
that contributes to the singularity in the k = (1, —1, 0) direction. As discussed earlier, such 
degeneracy gives rise to a logarithmic singularity in 7. 

A similar qualitative situation, for some different lattice theory [e.g. the FCC theory 
of (|7.3|) ], is shown in £g. Ipj. The same argument now goes through with the direction 
k = (1, 0, 0) and the pane of symmetry p^ = 0. This is the origin of the spikes in the FCC 
case of £g. Qc. 

To avoid spikes, look for lattice theories such that flp does not have lines of local maxima 
and minima in [p^, Py) that simultaneously he on the same plane of symmetry. To show that 
this is possible, I will give an example from massless scalar QED on a simple cubic lattice 
(where momenta are |pj| < tt). For scalar QED, the problem can be reduced to searching 
through the space of Dp that are hnite linear combinations of sin^(n-p/2) for integer vectors 
n. This is because one can generate any such D^, 

= I^CnSin^ , (7.10) 
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by a gauged version of the local scalar interaction 


i/scalar = E iHxl" + E E 

X X n ^ 

where If is conjugate to 0. So I will therefore hrst construct an f2p that works and then 
leave it as an exercise to the interested reader to construct the corresponding Hamiltonian. 
The space of linear combinations of sin^(n • p/2) is equivalent, by trigonometric identities, 
to the space of parity-even polynomials in cos pi and sinpj. 

A simple example of an flp that avoids spikes is 


= 4 - 2P(1 + cos Pi). 


(7.12) 


This example is a bit degenerate: the boundary of the Brillouin zone |pi| < vr is a constant 


energy surface, Hp = 4. A contour plot of Hp for hxed Pz is shown in £g. 


There 


is a maximum along planes of symmetry, but it is a degenerate maximum. The nominal 
conditions for a singularity in the damping rate are met but only at the boundary, and at 
the boundary the group velocity v corresponding to ([7.12| ) vanishes. The explicit factors of 
u* and in the formula (^.13| ) for 7 then compensate the divergence of the (5(v ■ k), and 
there are no singularities. 
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Here’s an example that’s degenerate in a different way: 

= {(1 — cos Pa;) {cospy — cosp^s)^! X {cyclic permntations} . (7-13) 

This example has been constructed to give fip = 0 along all the planes of symmetry, and 
a fixed pz slice is shown in fig. ^3. This example by itself doesn’t have the right p —> 0 
behavior for the continuum theory, but we can now obtain an interesting non-degenerate 
example by an appropriate superposition of the two degenerate ones. For example, a Pz slice 
of 

nl = + 5 isjpf (7.14) 

is plotted in fig. 


VIII. THE PLASMA FREQUENCY AND DEBYE MASS 


Since all the machinery was put in place in section |III A| , it’s interesting to look at some 
other quantities, characteristic of classical thermal gauge theories, that are not as directly 
related to the topological transition rate as the damping coefficient 7. In this section, I 
shall examine the plasma frequency and Debye screening length at leading order in weak 
coupling. 

The plasma frequency mpi is the frequency of small-amplitude, propagating, colored 
waves in the long wavelength limit.(In quantum theories, it is also known as the plasmon 
mass.) It can be extracted from the general hard-loop equation ( p.9|) by linearizing in 
and taking the k 0 limit for fixed oj. It’s convenient to work in a covariant gauge, where 
this limit gives 

" duQ 


= -2CAg^ 


dQ 


v’-y^A^. 


( 8 . 1 ) 


This concept certainly makes sense perturbatively. It is not clear (to me at least) whether there is a 
precise, useful, gauge-invariant, non-perturbative dehnition of the plasma frequency in non-Abelian gauge 
theory. 
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If the underlying theory has at least cubic symmetry, this can be rewritten as 




( 8 , 2 ) 


where 


2 2ri 2 f “"'0 I 12 

"^pl = -3<^A^ / — |v| . 

Jn 


dn^ 

/ p d^ 


(8.3) 


For a classical theory, this is 


ml = ICp^g^T f 

Jp 




(8.4) 


For the Kogut-Susskind Hamiltonian on a simple cubic lattice, the dispersion relationship 


( b-l|) gives 


ml, = 0.08606 Cas’ 


(8.5) 


in lattice units. 

The plasma frequency is in fact related to the physics of topological transitions. When 
the system crosses the potential energy barrier for such transitions, it oscillates many times 
back and forth across the barrier for each net transition [Q. These oscillations are small 
amplitude, and their oscillation frequency is the plasma frequency. For comparison, a typical 
SU(2) simulation by Ambjprn and Krasnitz IQ has couplings given by /^l = 4:/g‘^ = 12, for 
which (|8.5|) predicts a plasma oscillation period of 27r/mpi = 26.2. Their simulation has time 
steps of 0.05 in lattice units. 

Debye screening is the screening of static electric fields in a plasma. It can be seen from 


the small-amplitude, zero frequency behavior of (|3.9|) , which gives 


= mlA^ , 


( 8 , 6 ) 


where 


ml = -2CAg^ [ 


driQ 

dVt 


(8.7) 
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For a classical theory, 


ml = 2CAg‘^T [ 




( 8 . 8 ) 


For the Kogut-Susskind Hamiltonian on a simple cubic lattice, this gives 


ml = 0.50546 CaQ^ 


md = — 
TT^ 


(8.9) 


in lattice units. In fact, there is an analytic result 1110 for the integral, and (|8.9|) is more 
precisely 

o 4 


(18 + 12^2 - lOVs - 7VQ) [k(^{ 2 - VsfiVs - V2f)] Ca/ , (8 .10) 


where K{z) is the complete elliptic integral of the hrst kind. 

In the weak coupling limit, the inverse screening distance l/rud should be large compared 
to the lattice spacing. For comparison, the weakest coupling simulations in ref. have 
/3 l = 14, for which ( |0| ) gives l/m^ = 1.86 lattice spacings. So there is one physical quantity, 
at least, for which these simulations are only marginally in the weak-coupling limit. The 
implications for measuring the weak-coupling behavior of the topological transition rate are 
unclear. 


IX. CONCLUSION 

In this paper, I have presented a procedure for extracting the weak-coupling limit of 
the real topological transition rate from rates that could be measured in simulations. The 
concerns about rotational invariance originally raised by Bodeker et al. are crucial, and 
a precise measurement of the rate requires a search for an effectively rotational-invariant 
lattice Hamiltonian. Unlike the rotational improvements familiar to Euclidean studies, sim¬ 
ply improving the rotational invariance of higher and higher derivative terms in the small 


is I]/47r where E is given by eq. (A.5) of ref. [Q, except that equation has a typographic error: 
the elliptic function K should be squared 
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momentum limit is not enough: one must also introduce a UV cut-off distance A ^ that is 
large compared to the lattice spacing, as discussed in section B 

There remain several interesting problems. One is to explicitly construct the improved 
Hamiltonian of section Another is to survey other, simpler Hamiltonians (which would 
not require numerical extraction of the aA —0 limit) to hnd one with relatively small 
anisotropy as measured by It would also be interesting to hnd pure gauge theories 

that avoid logarithmic spikes, like the scalar QED theories discussed in sec. Ivra And, 
of course, there is the ongoing numerical problem of actually extracting lattice topological 
transition rates in the hrst place and verifying that they scale as at sufficiently small 
coupling. 

In this paper, I have focused on extracting the real topological rate from lattice simula¬ 
tions where the only degrees of freedom are standard gauge link variables. An interesting 
alternative suggested by Hu and Miiller [§] is to instead introduce additional particle degrees 
of freedom that are described by continuum positions x(t) and classical non-Abelian charges 
and which interact with the lattice helds. The purpose of these particles is to mimic the hard 
degrees of freedom of the real quantum theory and induce, in the lattice theory, the correct 
long-distance physics of the real theory. I would just like to point out that any proposal 
of this type will require very careful attention to orders of limits, because there are in fact 
two contributions to the effective long-distance theory: that of the new particles, and that 
of the short-distance modes of the lattice gauge helds. The latter gets large as the lattice 
spacing is made small and will be anisotropic. So two non-commuting limits are required 
of any such scheme: the limit of small lattice spacing, and a limit where the coupling to 
the additional particles is somehow made large so that their contribution to long-distance 
physics dominates. 

I am particularly indebted to Dam Son for explaining to me Blaizot and lancu’s formula¬ 
tion of hard thermal loops and for helping me understand how to generalize it. I also thank 
Steve Sharpe, Rajamani Narayan, Greg Moore, Larry McLerran, Misha Shaposhnikov, and 
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especially Larry Yaffe for a variety of useful conversations. This work was supported by the 
U.S. Department of Energy, grant DE-FG03-96ER40956. 


APPENDIX A: DIAGRAMMATIC DERIVATION OF 7*^ 


In this appendix, I will verify the formula ( |3.13|) for 7 *-^(k)—and in particular its inter¬ 
pretation of V as the group velocity—by applying diagrammatic methods in a wide class of 
theories. Specihcally, I will discuss the calculation of the imaginary part Im 11 of the gauge 
boson self-energy starting from one-loop Euclidean diagrams of the form of hg. ii, where 
the line in the loop represents hard excitations, be they gauge particles, scalars, or whatever. 
(Diagrams of the form of hg. will not give a contribution to Imll because they cannot 
be cut.) I will assume that tree-level propagators of the theory are of the form (pq — Dp)“^. 
Color indices, which are just responsible for the overall factor of Ca in the hnal result (|3.13|) 
for 7 , will be ignored. My hrst assertion is that the thermal contribution to Imfl*-^ is 

Imn*^'(o;,k) = (Ai) 


where u® is the distribution function for bosons; 


kl*i(o;,k) = ^ 


n. 


ab 


'p 


(a) ('1 4- 

2 ^ mi' - . 


2D 


(b) 

p+k 


(A 2 ) 


is the integrated, squared amplitude for absorption of a gauge boson by the thermal bath; 
and AI* = Al*(iL, P) is the vertex shown in hg. |l|. The sum above is over diherent species 
and polarizations a, b of the hard particles, and the capital letters K and P will henceforth 
stand for the “four-vectors” (a;,k) and (Dp,p). 1 ± Up represents the hnal state Bose or 
Fermi enhancement and should be replaced by just Up for a classical theory. 
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FIG. 5. Schematic contributions to the one-loop self-energy of soft gauge bosons. The loop denotes hard 
excitations of the theory. 


([Al|) is a simple consequence of the thermal optical theorem Rather than plodding 
through the one-loop case step by step here, I will give a quick pictorial argument for readers 
who are unfamiliar with it. By the contour trick for doing Euclidean hnite-temperature sums, 
or by simply starting directly with a real-time formulation of thermal perturbation theory, 
the thermal contribution to the self-energy diagram of £g. ^ is given by the diagrams of 
fig- P, which represent forward scattering of the soft excitation off of other excitations in the 
thermal bath. The imaginary part corresponds to cutting these diagrams and putting the cut 
line on shell. For the retarded self-energy, where uj is treated as a;-l-ie, this cutting yields the 
net decay rate (actually, the rate times —2a;) of excitations: —Im If is the difference A_ — 
shown in hg. ^ of the square amplitudes A- for destroying soft quanta and A^ for creating 
them. By the principle of detailed balance (which is indeed obeyed diagrammatically), these 
amplitudes must satisfy 


n^A+ = A_{1 +n ^), 


from which one then has 

ImB = -{A_-AA = . 

nuj 


(A3) 


(A4) 


37 



Phard 


Phard 


Phard 


Phard 




FIG. 6. A picture of the thermal part of fig. | 
The external hard lines are on shell. 


, as forward scattering off of particles in the thermal bath. 



FIG. 7. A picture of —Im 11 as corresponding to the difference between absorption and emission processes 
of soft quanta by scattering of hard particles in the plasma. 


So return to (^). Taking the limit that both uj and k are small and the limit that 
a; -C fc gives 


Imll*-^ ~ —27ij3uj E /mi±^p) 




< 5 (va ■ k) 


= 2^^ E ^ -7q 5-^K-k) 


(AS) 


The sum over flavors is now restricted to pairs of flavors that have degenerate dispersion 
relationships. The result (|3.13|) , or its generalization to other theories, will now follow if 
M.ah can be replaced by 5fV(r2p)5afe = 2gVL.pVa5ab'- that is, if soft magnetic gauge bosons 
indeed couple to group velocity. I’ll argue this by surveying a few instructive examples. 

As the first class of theories to check, start with a continuum scalar QED theory with a 
Lagrangian with scalar interactions of the form 
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^scalar ^ (dt^PaY{dt<Pa) “ /ab(D) (^k , 


(A6) 


where a and b are flavor indices. For the one flavor case, everything is trivial. The dispersion 
relationship is flp = /(p), and the trilinear conpling with the photon is g(j)*(j)A ■ V/; so Ai 
has the desired form. 

The case of mnltiple flavors is slightly less trivial and is worth nnderstanding becanse 
it is analogons to the case of mnltiple branches in lattice non-Abelian gange theory. Let 
U{p) be the nnitary matrix that diagonalizes /(p) in flavor space, so that F = fU is 
diagonal. It’s diagonal elements Fa are jnst the flp for the different branches. In this energy 
eigenbasis, the trilinear conpling is not in general diagonal: 

M = U^{Vf)U = VF + U^{VU)F + F{VU^)U = VF + [F^VF, F], (A7) 

so that 


Mab = yPaSab + {U^WU)ab{F, - Fa) . 


(A8) 


Bnt the degeneracy condition on the flavor snm in ( |A5|) eliminates the second term above, 
and so one obtains 


7*^(k) = -27r^/ 


driQ 

dVt 




(A9) 


This is indeed the trivial generalization of (|3.13|) to inclnde a snm over branches. 

Now consider scalar QED on a spatial lattice, with the gradient energy given by some 
snm of local terms of the form 


|0(x) - Fxy0(y)|% (AlO) 

where Fxy is the prodnct of U(l) link matrices along some path from x to y. (The canonical 
case wonld be that x and y are nearest neighbors and the path is the single link between 
them.) I am only interested in the conpling in pertnrbation theory and only in the 

limit of soft momentnm for A. In this limit, Fxy ~ 1 + 5 '(x — y) ■ A and the conpling is 
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minimal. That is, the coupling has the form ( |A^) in momentum space in the soft photon 
limit, and the previous discussion applies. 

Similar considerations hold for non-Abelian lattice gauge theories. Schematically, write 
A = Asoft + Ahard where Asoft are the soft-momentum modes and Ahard the hard momentum 
modes. Under soft gauge transformations G, we have Agoft —>■ — g~^G'VG~^ and 

Ahard —^ GAhardG”^- Gauge invariance under soft gauge transformations then implies that 
the AhardAsoftAhard coupling must be minimal in the soft limit. That is, it again has the 
form of ( |A^) in the soft limit with 0 representing Ahard and D = V — where T“ 

are the (anti-Hermitian) adjoint-representation generators. 

As a check, it is easy to explicitly verify that the soft-hard-hard coupling of the Kogut- 
Susskind Hamiltonian on a simple cubic lattice indeed has A4ab = g'^{^p)Sab, where a and 
b are the polarizations of the hard gluons. 


APPENDIX B: NUMERICAL INTEGRALS FOR try 

For the sake of completeness, I will explicitly write down the two-dimensional integrals 
corresponding to (|3.13|) that I did numerically to produce figs, ^b and c. 


1 . Simple cubic case 

vp 6{1 — I sinp^ 


/ dp^dpy Tt-f 

47r2|k| J-^ ^ Op 


COSPz 


where the Pz sum is over two values pf with 


sm Pz = —— sm Px — — sm py 
kz kz 


sm 


cospz = ±\Jl- sin Pz , 
^ (?) ^ l(l-cosp^). 


(Bl) 


(B2a) 

(B2b) 

(B2c) 


Op and V are given by (|5.1|) and (|5.2|) . 
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2. FCC case 


tr 7 (k) = 


87r2|k| 




r+27r 


/- 27 r 


dpxdpyJ2 


|v±| 


2 \/l 


r 2 ± Pz 


B cos — A sin 

xe{A^ + -c^)e (i 


, Pz 

smi 2 


(B3) 


where the f2 sum is over the two dispersion relations (|7.7|) ; Cp and v± are given by (|7.5|) and 
( [7.8|) with f2p = f2±; the pz sum is over two values pf with 


and 


BC ± A^A^ + B^-C^ 
A^ + B^ 


'Pz 

sm 1 y I = 

'Pz\ AC T BVA^ + B^-C^ 

' y j “ ^42 + 52 


(B4a) 

(B4b) 


A = kz 


B = kz 


sin 1 + sin 

cos I I) + cos (I 


> • iPx\ fPy\ f . (Py\ (Pa 

C = —kr sm — cos — — fc,, sm — cos 


(B5) 

(B6) 

(B7) 


The ±(=f) sign in (p^) is unrelated to that in (|7.7|) 
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